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Abstract. We study differential operators on an elliptic curve of order higher than 2 which 
are algebraically integrable (i.e., finite gap). We discuss classification of such operators 
of order 3 with one pole, discovering exotic operators on special elliptic curves defined 
over Q which do not deform to generic elliptic curves. We also study algebraically integrable 
operators of higher order with several poles and with symmetries, and (conjecturally) relate 
1^ ' them to crystallographic elliptic Calogero-Moser systems (which is a generalization of the 

d , results of Airault, McKean, and Moser). 
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O ! To Igor Moiseevich Krichever on his 60th birthday 

^ ! 1 Introduction 

In this paper we study differential operators 

L = 9" + 02(z)a"-2 + ... + a„(^) 

^ \ which are algebraically integrable (i.e., there exists a nonzero differential operator M of order 

^ ' relatively prime to L such that [L, M] = 0). Such operators were first studied in [2] and became 

a focus of attention since the seventies, as they provide explicit solutions to the Gel'fand-Dickey 
hierarchy (in particular, the KdV hierarchy for n = 2 and the Boussinesq hierarchy for n = 3; 
see [H [7] and references therein) . A general classification of such operators was obtained in [13] . 

We are interested in making this classification more explicit in the special case when the 
coefficients ai{z) are meromorphic functions on an elliptic curve E. For instance, in the simplest 
nontrivial case n = 2 and a single pole, it is well known that the only algebraically integrable 
operator, up to equivalence, is the Lame operator 

L = — m{m + l)p{z), 

where m is an integer; its algebraic integrability was discovered by Hermite. In the case n = 2 
and multiple poles, the answer is much more interesting (see [H Section 4.1], as well as fiQ\ [TT] 
and references therein). 



*This paper is a contribution to the Special Issue "Relationship of Ortliogonal Polynomials and Spe- 
cial Functions with Quantum Groups and Integrable Systems". The full collection is available at 
|http://www.emis.de/journals/SIGMA/QPSF.html| 
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We study the problem of classification of algebraically integrable operators L for n > 2. It 
turns out that already in the case n = 3 and a single pole, the situation is much richer and more 
complicated than for n = 2; in particular, there exist algebraically integrable operators L of third 
order with one pole defined on an infinite family of special elliptic curves over Q, which do not 
deform to operators (with one pole) on a generic elliptic curve. We provide a list of third order 
algebraically integrable operators with one pole which is conjecturally complete, and state some 
results and conjectures concerning operators with several poles. In particular, we conjecture that 
in the special case of operators with symmetries, algebraically integrable operators are described 
in terms of the classical crystallographic elliptic Calogero-Moser systems introduced in [8]. 

The paper is organized as follows. In Section [2l we give an exposition of the general theory 
of algebraically integrable operators, in particular those on an elliptic curve; the results here 
are mostly well known, but for reader's convenience we give an exposition based on differential 
Galois theory similar to one in [1]. In this section we also propose a general conjecture on the 
classification of operators with one pole. In Section [3l we present computational results for 
third order algebraically integrable operators with one pole, and give a conjectural classification 
of such operators. Finally, in Section H] we discuss operators with several poles, and state 
a conjecture on the connection with the systems of [8]. 

2 The general theory of algebraically integrable operators 

In this subsection we review the basics on algebraically integrable ordinary differential operators. 
Most of this material is well known; we refer the reader to [131 IISl IISl EZl ttOl HI 13 E] and 
references therein. 

2.1 Definition of algebraic integrability and meromorphicity of coefficients 

Consider the differential operator 



where ai{z) are smooth functions on some interval in M. Note that the coefficient ai can be 



Recall that L is called algebraically integrable (or algebro-geometric) if there exists a nonzero 
differential operator M of order relatively prime to L such that [L, M] = [2]. Note that up to 
scaling M is necessarily monic (i.e., has leading coefficient 1). 

For example, if n = 2 then L = d"^ + u(z), and if L is algebraically integrable then u is called 
a finite- gap potential. 

Theorem 2.1 ( [16\ Theorem 6.10], see also [I^). If L is algebraically integrable, then ai{z) 
extend to meromorphic functions on the complex plane. Moreover, the order of each pole of the 
function Oj in C is at most i, for i = 2, . . . ,n; in other words, the operator L has regular {or 
Fuchsian) singularities in C. 

2.2 The indices of L 

Let ai{z) = hiZ~^{l + 0{z)) near z = 0. Then by rescaling z the operator L can be degenerated 
into the operator with rational coefficients 



L = 9" + ai(z)a"-^ + --- + a„(z) 




& 



+ 622"^5"-2 + • • • + bnZ-'^ = 0. 
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Consider the differential equation 
Lqz"^ = 0. 

This equation is equivalent to the algebraic equation 

PL{m) = 0, 
where 

PL{m) = m{m — 1) • • • (m — n + 1) + b2m{m — 1) • • • (m — ?i + 3) + • • • + 6„. 

Let rrij, j = 0, . . . ,n — 1, be the roots (with multiplicities) of the polynomial Pl, i.e. 

-Pl(s) = (s - "lo) • • • (s - 

The numbers nij are called the indices of L at 0. They are arbitrary numbers satisfying the 
relation 

n-l 

rrij = n{n — l)/2. 

i=o 

Obviously, the indices uniquely determine the coefficients bi. 
Similarly, one defines the indices of L at any point zq. 

Example 2.2. The indices of L at a regular point are 0, 1, . . . , n — 1. 

2.3 Algebraic integrability of homogeneous rational operators 

Proposition 2.3. The operator Lq given by formula (j2.ip is algebraically integrable if and only 
if the indices mj are integers which are distinct modulo n. 

Proof. It is easy to see (see e.g. [3]) that the operator Lq is algebraically integrable if and only 
if the equation Lqi^ = fi'^ip admits a Baker- Akhiezer solution of the form F{jjLz), where 

F{x) = e^Q{l/x), 

and Q is a polynomial such that (5(0) = 10 Solving the equation by the power series method, 
we see that this happens if and only if mj are distinct modulo n (and thus represent each residue 
class exactly once). ■ 

Corollary 2.4. If L is algebraically integrable then the indices mj{zo) of L at every point zq 
are integers which are distinct modulo n. 

Proof. Suppose L is algebraically integrable. Since by Theorem 12.11 the commuting opera- 
tor M has regular singularities at zq, its rational degeneration Mq is well defined, and commutes 
with Lq, so Lq is algebraically integrable. So the result follows from Proposition 12.31 ■ 

When rrij are integers, we will order them in the increasing order, tuq < mi < ■ ■ ■ < m„,_i. 
It is also convenient to introduce the "gaps" qj := mj — nij-i, j = 1, . . . , n — 1, which clearly 
determine mj. In the integrable case, these are nonnegative integers not divisible by n. 

^E.g., if Lo is algebraically integrable, then it admits a homogeneous commuting operator Mo of relatively 
prime order m such that = AIq, and the system of differential equations LqV' ~ f^"i^, Motp — can be 

reduced to a first order scalar equation, which has a solution of the required form by Euler's formula. 
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2.4 Algebraic integrability of differential operators on elliptic curves 

Let r C C be a lattice, and E = C/T he the corresponding elliptic curve. Assume that ai{z) 
are rational functions on E (i.e., elliptic functions). 

Theorem 2.5. The following conditions on L are equivalent. 

(i) L is algebraically integrahle. 
(ii) The monodromy of the equation 

LiIj = Xijj (2.2) 

around every pole of L in E is trivial for any eigenvalue A € C 

[Hi) The monodromy group of equation ()2.2p is upper triangular in some basis. 

(iv) For generic A {i.e., outside of finitely many values), equation (j2.2p has a basis of solutions 
of the form 

where 9 is the first Jacobi theta- function. 

Proof, {i) =^ (ii). Since L is algebraically integrable, by Corollary 12. 4| the indices of L at 
every pole are integers. Hence the monodromy matrices of equation (12. 2|) around the poles of L 
are unipotent. Also, it follows from [3] that the differential Galois group of ()2.2p for generic A 
is an algebraic torus. Since monodromy matrices belong to the differential Galois group, and 
since every unipotent element of a torus is trivial, we conclude that the monodromy matrices 
of (12. 2p around the poles of L are trivial for generic, hence for all A. 

(ii) =^ (Hi). If (ii) holds, the monodromy group of ()2.2p is Abelian, hence is upper 
triangular in some basis. 

(Hi) =^ (if). Assume (Hi) holds. Since L has regular singularities, the differential Galois 
group of equation (12. 2p is the Zariski closure of the monodromy group. Hence the differential 
Galois group of equation (12. 2p is triangular as well. But it is shown in [3] that the differential 
Galois group of (j2.2p is reductive for generic A. Hence this group is Abelian, and is a torus for 
generic A. Thus there is a fundamental system of solutions of the form ()2.3p (see [3]). 

(iv) =^ (i). If (iv) holds, then the monodromy of equation (j2.2p around poles is trivial for 
generic A. Hence it is trivial for all A, and the monodromy group, hence the differential Galois 
group of ()2.2p is Abelian. So by [3], there exists a nonzero differential operator M of order 
coprime to the order of L such that [L, M] =0. ■ 

Remark 2.6. 

1. This theorem is similar to Theorem 5.9 in (5], which goes back to [4]. 

2. A similar theorem, with the same proof, holds in the trigonometric (nodal) and rational 
(cuspidal) case, i.e., when the coefficients of L are rational functions on the nodal or cuspidal 
curve of arithmetic genus 1 which are regular at infinity. More precisely, in the rational case, 
since the singularity at infinity is irregular, we must add the Stokes matrix at infinity to the 
monodromy group, and we should also replace 6{z) with sin(z) and z in the trigonometric and 
rational cases, respectively. We note that in the trigonometric and rational case, the implication 
(iv) =^ (i) of Theorem 12.51 was proved in |19j . 

2.5 Operators with one pole 

In this subsection we will consider the special case when L has only one pole, at the point £ E. 
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2.5.1 Second order operators with one pole 

Let n = 2, and let L have a unique pole at (the simplest nontrivial case). In this case, up to 
an additive constant, the operator L has the form 

L = d"^ + ap{z), 

where p is the Weierstrass function of E (the Lame operator). Local analysis near (i.e., the 
condition that the local monodromy is trivial) shows that algebraic integrability of such L implies 
that a = —m{m + 1), where m is a nonnegative integer. Conversely, it was shown by Hermite 
that if a = —m{m + 1) then L is indeed algebraically integrable. Namely, the triviality of the 
monodromy of (|2.2|) near is easy to see by noting that the operator L is Z2-invariant, while 
the indices at are —m and m + 1, whose difference is odd; thus, the algebraic integrability 
of L follows from Theorem 12.51 

2.5.2 The algebraic integrability locus 

Assume now that we have fixed the indices mj distinct modulo n (and thus the coefficients 6j). 
Then the set of possible operators L is parametrized by the Laurent coefficients Oj^ of ai{z) of 
nonpositive degrees —k, k < i, and also by the coefficients g2, 93 of the differential equation 

for the Weierstrass function p of the elliptic curve E. Note that we have the C*-action rescaling 
the lattice F, with respect to which these parameters have the following degrees (or weights): 

deg(52) = 4, degigs) = 6, deg(aifc) = i - k. 

Thus, we will think of these parameters as homogeneous coordinates on the weighted projective 
space with these weights, and define the algebraic integrability locus (for each choice of 
indices m) as a subset of this weighted projective space. 

2.5.3 The cyclically symmetric operators 

For every choice of the indices m, there exists a unique operator L such that all the coefficients 
Oj/c equal zero. Let us denote this operator by L{0). 

Proposition 2.7. For n > 2, the operator L{0) is algebraically integrable {for any indices 
distinct modulo n) in the following cases: 

(i) n = 2 [the Lame operator any elliptic curve); 
[a) n = 3, 6, (72 = (equianharmonic elliptic curve); 
(Hi) n = 4, (73 = [lemniscatic elliptic curve). 

Proof. In these cases, the operator L(0) has a symmetry under the groups Z3, Zg, and Z4, 
respectively. This symmetry easily implies the triviality of monodromy at 0. ■ 

Remark 2.8. Case (i) is well known and due to Hermite. Case {ii) was done by Halphen 
[m p. 571] in the case of equal gaps; his proof easily extends to the case of general gaps and 
general n. For multivariable generalizations of Proposition YI7j\ see [8]. 

Note that the operators L(0) — A, where iv(0) is as in Proposition 12.71 {i)~{iii) are the only 
operators L which are symmetric under Z„ (where n is the order of L), acting by 
We will call such operators fully cyclically symmetric (or just cyclically symmetric if there is no 
ambiguity (e.g., when n is a prime)). 
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Conjecture 2.9. For any n>2, there exists N G Z+ such that if qj > N for all j = 1, . . . , n — 1, 

then the only algebraically integrahle operators L with gaps qj are fully cyclically symmetric. 

Remark 2.10. 

1. Coniecture l2.9| in particular, claims that algebraic integrability for large enough gaps takes 
place only for n = 2, 3, 4, 6. 

2. As we have explained, for n = 2 the Conjecture 12.91 holds with = 0. 

3. Conjecture 12.91 is open even for n = 3. It is supported by computational evidence and 
partial results described in the next section. 

A similar conjecture can be made in the rational case. Namely, assume that aj(z) G C[z~"'^]. 
In this case, we have fully cyclically symmetric operators Lq — A for any n > 2, which are 
algebraically integrable. 

Conjecture 2.11. Conjecture 12.91 holds in the rational case. 

Remark 2.12. In the case n = 2, the differential equation Lif) = Xip is confluent hypergeomet- 
ric, and Conjecture 12.111 is well known to be true (with A^ = 0). 

2.6 The classification of elliptic finite gap potentials 

In this subsection we give a proof of the classification theorem of finite-gap potentials on elliptic 
curves with arbitrary number of poles |im Bl[TT]: this proof is based on differential Galois theory 
and follows [H Section 4.1]. 

Theorem 2.13 ([H]). Let u{z) be a rational function on an elliptic curve E, which is a finite-gap 
potential {i.e., the operator L = d"^ -\-u is algebraically integrahle). Then there exist nonnegative 
integers mi, . . . ,mN and points zi, . . . , zn & E satisfying the equations 

'^mj{mj + l)p^^''~'^^Zi- Zj) = 0, i = l,...,N, s = l,...,mj, 
such that 

N 

u{z) = — rni{mi + l)p(z — Zi) + const. 
1=1 

Moreover, any potential of this form is finite-gap. 

For the proof of this theorem, we will need the following classical lemma from the theory of 
differential equations. 

Lemma 2.14. Let m G Z+, and 

oo 

u{z) = -m{m + l)z~^ + ^ CjZ^ G 
Then the equation 

admits two linearly independent solutions in C((z)) for all X & C if and only if ci = c^ = ■ ■ ■ = 
C2m-1 = 0. 
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Proof. We need to show that the given condition is equivalent to the existence of a solution 
of the form X^„>o ^n-^""*" with ag = 1. If all odd-numbered coefficients Cj are zero (the Z2- 
symmetric case), then the required solution clearly exists. Otherwise, the obstruction to the 
existence of such solution is a certain homogeneous polynomial P of Cj and A of degree 2m + 1 
(where deg(cj) = j + 2, deg(A) = 2). Suppose that s is the smallest integer such that C2s-i 
is nonzero. Then it is easy to compute that the polynomial P has degree m — s with respect 
to A, and its leading term is a multiple of C2s-iA"*~*. Thus, this polynomial is nonzero, and the 
required solution does not exist. On the other hand, if s > m + 1, then P = 0, and the required 
solution exists. ■ 

Proof of Theorem 12.131 Theorem 12.11 implies that the poles of u must be exactly second 
order (it is clear that if u has a first order pole then ()2.2p does not admit a basis of meromorphic 
solutions). Moreover, it follows from Corollary 12.41 that the singular part of u at each pole Zi, 
i = 1, . . . , N , is —mi{mi + l)(z — Zi)~'^, where mi is a nonnegative integer. Thus, we find that 
in the algebraically integrable case 

N 



u{z) = — rni{mi + l)p{z — Zi) + const. 



i=l 



Let us now show that the conditions for algebraic integrability in terms of rrii , Zi are exactly 
as stated in the theorem. According to Theorem 12.51 the condition for algebraic integrability 
is that the monodromy of (j2.2p around each pole of u is trivial. So the theorem follows from 
Lemma 12.141 ■ 

In particular, if all rrii = 1, we obtain the following well known result: 

Corollary 2.15 (|lj). The potential u = —2^ p{z — zi) is algebraically integrable if and only 
if (zi, . . . , zn) is a critical point of the elliptic Calogero-Moser potential 

U{zi,...,zn) = ^ p{zi-Zj). 
l<j^i<N 

Remark 2.16. The same method can be used to rederive the classification from |10j of trigono- 
metric and rational finite-gap potentials which are bounded at infinity, i.e., potentials on the 
nodal and the cuspidal curve of arithmetic genus l| It leads to the same answer, with p{z) 
replaced by and respectively. 

Note that (as explained in [10] and references therein), in the rational case we have the 
identity 

mj(mj -|- 1) = m{m + 1) (2-4) 

i 

for some m G Z+. This identity comes from the fact that a rational potential 

can be degenerated into 

Y,imi{mi + 1) 



Uo{z) 



Z2 



^As explained in [TT], boundedness at infinity is in fact automatic. 
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by rescaling z (so that algebraic integrability of 9^ + u implies algebraic integrability of 9^ + uo)- 
We see that, surprisingly, identity (|2.4|) follows automatically from the existence of solutions of 
the equations 

3 Third order operators with one pole 

3.1 The general setup 

Let us now consider in detail the case n = 3 with one pole at 0. Up to an additive constant, the 
operator L in this case looks like 

L = + (ap{z) + c)d + + ep(z)), (3.1) 

where a,b,c,e S C. 

As explained above, a necessary condition for algebraic integrability is that the indices niQ < 
mi < m2 at are integers pairwise distinct modulo 3. So we have mo + mi + m2 = 3, and 

a = m^mi + mom2 + mim2 — 2, b = -mim2m^. 

As above, it is convenient to introduce the "gaps" q = qi = mi — mo, r = q2 = m2 — mi (then 
mo = 1 - (2g + r)/3, mi = I + {q - r)/3, ms = 1 + (2r + q)/3). 

Our goal is to determine the algebraic integrability locus for each set mo, mi, m2 in terms 
of the homogeneous coordinates c, e, g2, gs- 

Recall that the j-invariant of E is defined by the formula 

1728^^ 
J (^) - ' 



9i - 27gi 

In particular, the equianharmonic (i.e., Za-symmetric) elliptic curve E has 52 = and j = 0. 
Thus, by Theorem 12.71 the point c = e = 52 = belongs to the algebraic integrability locus for 
any q, r. 

Note that the parameters o, h of the operator — L*, where L* is the adjoint operator to L, 
are given by the formulas a' = a, 6' = — 6 + a (and also c' = c, e' = — e). So the indices of L* 
are 2 — mj, i = 0, 1, 2, and hence the gaps q and r are interchanged under passing to the adjoint 
operator. On the other hand, it is clear that L* is algebraically integrable if and only if so is L. 
So it suffices to consider the case q > r. 

Also note that the gaps cannot be divisible by 3, and must be equal modulo 3. So we can 
write q = r + 3k, where /c > is an integer. 

We have a basis of solutions of equation (|2.2p of the form 

V^,(z) = z"Hl + o(l)), ^^0. 

Obviously, ^2 is single- valued near 0; algebraic integrability of L is equivalent to the condition 
that ipo{z), i^i^z) are single- valued near 0, i.e., do not contain log factors. However, it turns out 
that even a weaker condition suffices. Namely, we have the following important proposition. 

Proposition 3.1. Iftl^i is single-valued near for all A, then so is 'i/'o? CLiT-d thus L is algebraically 
integrable. 

The proof of Proposition 13. II is based on the following well known lemma from linear algebra. 
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Lemma 3.2. // A, B are two square matrices such that AB — BA has rank at most 1 then 
A, B are simultaneously upper triangular in some basis. 

Proof. Without loss of generality, we can assume ker^ ^ (by replacing A with ^ — A if 
needed) and that A ^ 0. It suffices to show that there exists a proper nonzero subspace 
invariant under A, B; then the statement will follow by induction in dimension. 

Let C = [A, B] and suppose rank C = 1 (since the case rank C = is trivial). If ker A C ker C, 
then ker j4 is i?-invariant: if Av = then ABv = BAv + Cv = 0. Thus ker^ is the required 
subspace. If ker A ^ kerC, then there exists a vector v such that Av = but Cv ^ 0. So 
ABv = Cv ^ 0. Thus ImC C Imvl. So Im^ is S-invariant: BAv = ABv + Cv € ImA. So 
Im A is the required subspace. 

This proves the lemma. ■ 

Corollary 3.3. If A, B are two invertible square matrices such that ABA^^B^^ — 1 has rank 
at most 1 then A, B are simultaneously upper triangular in some basis. 

Proof. Let ABA~^B~^ -l = T, where T has rank at most 1. Then AB - BA = TBA, and 
TBA has rank at most 1. So the result follows from the lemma. ■ 

Proof of Proposition 13.11 Let A and B be the monodromy matrices of ()2.2p (for some base 
point) along the global cycles of the elliptic curve. Then ABA~^B~^ = U, where U is unipotent. 
Since ipi is single valued, U — 1 has rank 1. So it follows from Corollary [33] that the monodromy 
of (j2.2p is upper triangular, so the result follows from Theorem 12.51 ■ 

Using Proposition 13.11 and CAS "Maple" , we computed the algebraic integrability locus for 
small values of r and any value of q. In each case, we applied L — A to a generic function 

^-1 (1 + + f^,2 + ... + f^^.r ^ O (/•+!)), 

solved for /i, . . . , fr-i, and look at the coefficient of in the result. This is a polynomial 

in A, and by Proposition 13.11 the monodromy around is trivial iff the coefficients of this 
polynomial vanish; this gives equations in c, e, §2, §3. 

Remark 3.4. As explained above, homothety of the elliptic curve rescales c, e, g2, gs and A, 
giving a natural notion of degree: 

deg(e) = 1, deg(c) = 2, deg{g2) = 4, degigs) = 6, deg(A) = 3, 

and the constraint polynomial is homogeneous of degree r. 

The results of our computations are presented in the next subsection. 

3.2 Results of computer calculations 
3.2.1 r = 1 

The coefficient of in the image of L — A is e, and thus the operator is algebraically 

integrable iff e = 0; this gives a 1-parameter family of integrable operators on each elliptic 
curve. The case g' = 1 is trivial (L = 5^ + cd), while the case q = A was considered by Picard in 
1881. 



10 



p. Etingof and E. Rains 



3.2.2 r = 2 

The coefficient of z™^"'"^, after solving for /i, is independent of A, and thus gives an equation 
relating c and e: 



3e' 



2 



Each curve thus has a 1-parameter family of integrable operators of this form. 

In particular, setting q = r, e = we get Example 1 of [18]. This is the third order operator 
commuting with the Lame operator 5^ — 2p. 



3.2.3 r = 4 

Now the coefficient is linear in A. The leading coefficient is a positive multiple of e, so e = 0; 
then the constant coefficient relates g2 and c: 

c = -^92. 

Thus each elliptic curve admits a pair of integrable operators of this form, except that when 
g2 = 0, the two operators coincide (with the corresponding cyclically symmetric operator). 
In particular, setting g = r, e = we get Example 2 of |18j . 



3.2.4 r = 5 

Again the constraining polynomial is linear in A. The leading coefficient implies 

3(7g2 + 35g + 46) 2 
"^""16(9 + 1)2(^ + 4)2'' 

at which point the constant coefficient factors, so that either c = e = or 

_ 27(4g2 + 20g + 25) 4 
64((? + l)4(g + 4)4 • 

Each elliptic curve thus admits four integrable operators of this form (or two modulo the sym- 
metry z — >• — z) except the equianharmonic case 52 = 0, where all these operators coincide with 
the cyclically symmetric operator. 

In particular, setting q = r, e = we get Example 3 of [18]. 



3.2.5 r = 7 

The constraint polynomial has degree 2, with leading coefficient 

5(g + 2)(g + 5) 
lU{q + l){q + 3)iq + 4){q + 6f' 

so that e must be 0. But this also eliminates the constant term, leaving only an equation 
relating c and §2- 

2_ 25(g + 2)2(g + 5)2 
12(2g + 7)2 

So we get two operators on each elliptic curve which coincide in the Zs-symmetric case 52 = 
(with the cyclically symmetric operator). 

In particular, setting g = r, e = we get Example 4 of [18] . 
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3.2.6 r = 8 

The constraint polynomial again has degree 2, and the leading coefficient implies 

3(191g'^ + 3056^3 + 17598^2 + 42992g + 38384) ^ 
" ~ 686(g + 1)2(^ + 4)2(^ + 7)2 ^ • 

Unlike in the previous case, however, one still has two conditions remaining, of weighted degrees 5 
and 8 respectively. The degree 5 condition naturally has a factor of e, but setting e = 
makes c = so that the remaining constraint is a multiple of §2', the coefficient is negative, so 
this makes 52 = 0. In the remaining case, the degree 5 condition gives a formula for §2'- 



92 



„den 
92 



where 



and 



gf^ ■= 155383^^ + 4972256^^ + 68978821^^ + 541706360g^ + 2632855228?" 
+ 8104425920?^ + 15416669104?^ + 16555419008? + 767835508 



^den ^ 470596(g + l)4(g + 4)4(g + 7)"(l9?2 + 152? + 277) , 



and the degree 8 condition can then be solved for gs/e^. For each ? = 8 + 3/c, A; > 0, one thus has 
a single j-invariant other than j = for which there exists an algebraically integrable operator: 



J 



-6912- 



;dcn 



where 



wit 



iH 



P8(g)^(lV + 152? + 277), 



ps{q) = 155383?^ + 4972256?^^ + 68978821?^ + 541706360?^ + 2632855228?" 
+ 8104425920?^ + 15416669104?^ + 16555419008? + 7678355008 



and 



•den 



(? + 7)(? + 6)(? + 2)(? + 1)(67?2 + 533? + 898) (67?^ + 539? + 922) 
(37?^ + 399?2 + 1344? + 1468) (37?^ + 489?^ + 2064? + 2692) 
(367?^ + 5115?^ + 23376? + 34828)(367?^ + 3693?^ + 12000? + 12724) 
(829?" + 14194?^ + 89097?^ + 242068? + 239236) 
(829?" + 12334?^ + 66777?^ + 156028? + 133156) . 



Note that j is finite for any ? since the factors of j'^'^'^ have positive coefficients. 

In particular, for e = the only solution is §2 = 0, which is shown for ? = r in Example 5 

of 



^We note that ps is the numerator of 32/6''. 
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3.2.7 r = 10 

The constraint polynomial is cubic in A, and the leading coefficient implies that e = 0; setting 
e = makes the polynomial even in A, so one has two remaining constraints, of degrees 4 and 10 
respectively. The degree 4 constraint can be solved for 52- 

_ 3(2069g^ + 41380g3 + 301017^2 _^ 941170^ + 1071464) 2 
~ 4400(g + 2)2(g + 5)2(g + 8)2 ' 

at which point the degree 10 constraint is (? times an equation for 33. Thus either e = c = g'2 = 
or 



53 
where 



„num 

.53 

„den ' 
53 



■= 96577g^ + 2897310q^ + 35259207^^ + 222299140g^ 
+ 764656215g2 + 1360455150g + 978817201 

and 

^den ._ 422400(g + 2f{q + 5f{q + Sf. 
Thus other than j = 0, the only possible j invariant is 

„-num 

= -995328^, 

where 

^•num ._ (2069^'' + 41380g^ + 301017g2 + 941170g + 1071464)^ 

and 

f-^"^ = (5g + 19) (5(7 + 31)(13g + 47)(13g + 83) (17,? + 73)(17g + 97)(19g + 59) 

X {19q + 131) (11^2 _^ iiQq _^ 239) (23^2 _^ 200g + 317) (23^2 _^ 260g + 617) . 



3.2.8 r = 11 

The constraint polynomial is cubic in A, and other than c = e = g2 = 0, there is no solution to 
the resulting four equations; one can solve the first three for c, g2, gs in terms of e, and plug in 
to the fourth equation, obtaining e^^ times a rational function which is negative for q > 8. Thus 
c = e = gf2 = is the only solution. 



3.2.9 r = 13 

The constraint polynomial is quartic in A, but again the first equation is e = 0, and eliminates 
half of the remaining equations. One thus has two additional equations which can be solved to 
give either c = 32 = or 52 /c2, 53 /c^ equal to specific rational functions of q. There is thus 
again a single surviving j invariant, which is given by the formula 

„-num 

^■ = -124416^, 
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where 

j-num ^ (g7^2 _^ gy^^ _^ 2014)2 (24727^^ + 964353g^ + 15225009^^ + 124224139g3 
+ 551142996^2 + 1258400208^ + 1155995968)^ 

and 

^•dcn ^ ^^3^ _^ 68)(13g + 101)(l9g2 + 265,? + 796) (19^^ _^ 229g + 562) 

X (83g2 + 1094(? + 2936) {83q'^ + 1064g + 2741) (47g^ + 924^^ + 5481g + 9532) 
X (47g^ + 909g2 + 5286g + 8824) (547g^ + 14649^2 ^ 1277589 + 360056) 
X (547g^ + 6684g2 + 24213g + 26876) (llf?^ + l43g + 332) . 

3.2.10 r = 14 

If e = 0, then c = (72 = 0; otherwise, one can solve the first three equations for c/e^, g2/e^, 

at which point the fourth equation is again a negative multiple of e^^, so there is no additional 

solution. Similar arguments apply to r = 17, 20. 

3.2.11 r = 16 

The first equation makes e = 0, so that one has three additional equations. If c = 0, then the 
next equation makes g2 = 0; otherwise, one can solve the first two equations for §2/0^, 93 /c^, at 
which point the remaining equation is a negative (for q > 13) multiple of c^, so no other solution 
exists. Similar arguments apply to r = 19, 22. 

Remark 3.5. Observe that the polynomials j'^^^ in the cases r = 8, 10, 13 split into many 
irreducible factors over Q, whose leading coefficients are either 1 or primes. Moreover, the 
constant coefficients of the factors for r = 10 are also primes, while for r = 8, 13 they are primes 
times a small (at most third) power of 2. The number-theoretic roots of this peculiar behavior 
are mysterious to us. We burden the reader with the unwieldy expressions of the j-invariants in 
the hope that someone would help us demystify it. 

Remark 3.6. The operators with q = r for r = 10 and r = 13 exist only for special values of j, 
so they are not present in [18], which deals with the case of generic j only. 

3.3 The classification conjecture 

On the basis of this data we make the following conjecture. 

Conjecture 3.7. For q > r > 14, there are no algebraically integrable operators L apart from 
the one which is Z,^- symmetric {i.e. c = e = 52 = 0). Thus, all the algebraically integrable third 
order operators L with one pole are the ones described in this subsection. 

This is a more precise version of Conjecture 12.91 for n = 3, claiming that in this case one may 
take N = 14. 

Here is a partial result in the direction of this conjecture. 
Proposition 3.8. 

(i) If r = 2>s + 1, where s > is an integer, then for any algebraically integrable operator L of 
the form (|3.ip . one has e = 0. 
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(ii) Assume that r = 3s + 1, where s > 1 is an integer. If L is an algebraically integrable 
operator of the form (|3.ip then g2 = if and only if c = {in which case L is cyclically 
symmetric). 

Proof, (i) We have already checked the case 5 = directly, so we may assume that s > 0. 

Algebraic integrability of the operator (j3.ip is equivalent to the existence of three linearly 
independent solutions of the differential equation 

d^f + (t^apitz) + t^c)df + {thp'itz) + t^ep{tz))f = f 

in C{{z)) for generic t; this is just the image of the original eigenvalue equation under a homothety 
z I— 7- t2 of scale t = A~^/'^. For t = 0, this has three independent solutions, each of which must 
deform to a solution for general t with the same asymptotics at z = 0. Let tpi{z) be the middle 
solution, 



^^^(l-r/3)fc(l + g/3)fcfc!- 

Then we need to be able to deform this to a solution of the form ipi{z) + tGi{z) + 0{t^), and 
we claim that this implies e = 0. Plugging 



Gi(z)=^c(fe)z-- 



i+3fc+l 

:.ytv 

k>0 



into the equation gives 

^[(3A; + l){3k + l-r){3k + q + l)c{k) - c{k - l)]z"'i+3k-2 



k>0 

^rrii+3fc— 2 

= - E 33^(1 -./3),(l + ,/3U. (^^ + 1)(3^ + 1 - r)i3k + , + l)c(.) - dk - 1) 



33fc(l-r/3)fc(l + (//3)fcA:!- 

Setting c'{k) = c{k)3^''{l - r/3)fc(l + q/3)kk\, we obtain 

{3k + l)(r -3k- l){q + 3k + l)c{k) - 3k{r - 3k){3k + q)c{k - 1) = e. 

It follows that c'{k) is a positive multiple of e for < A; < (r — l)/3. Since the equation for 
A; = (r — l)/3 reads —{r— l){q + r — l)c'((r — 4)/3) = e, this gives a contradiction unless e = 0. 

(ii) We have already computed the case s = 1 directly, so we may assume s > 1. Also, it 
follows from part (i) that e = 0. The perturbed solution, if it exists, will thus have the form 

iPi (z) + t^Gi {z) + t^G2 {z) + O [t^] . 

(Note that the odd degree terms vanish by symmetry.) To the same order, the differential 
equation reads 

Writing 

(z/3)3'^ 



{l-r/3)k{l + q/3)kkV 
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G2{z) = y d'{k)- — — — — 

^ ' Y (l-^/3)fc(l + g/3)fcA;! 

and substituting in, we find from the term that 

{2,k + 2)(r - 3A: - 2){q + 3/c + 2)c{k) - 3k{r - 3k){q + 3k)c{k - 1) 
= {{q-r)/3 + 3k + l)c 

and thus c'{k) is a positive multiple of c for < A; < (r — 4)/3. We also find that 

{3k + l)(r - 3k -l){q + 3k + l)d'{k) - 3k{r - 3k){q + 3k)d'{k - 1) 

= -k{r - 3k){q + 3k){q-r + 9k)c{k - l)c + k{r - 3k){q + 3k)C{k, q, r)g2 

where C{k,q,r) is positive when 4<3k + l<r<q, since the appropriate linear change 
of variables gives a polynomial with positive coefficients. We find by induction that d'{k) is 
a nonnegative linear combination of — and 52 for < A: < (r — 4)/3, while the equation for 
k = (r — l)/3 also tells us that d'{{r — 4)/3) is a negative linear combination of — and g2. 
Subtracting the expressions for d'{(r — 4)/3) gives a positive linear combination of — and g2 
which vanishes. In particular, if one of c and g2 is 0, so is the other. ■ 

3.4 The nodal and cuspidal cases 

The results of the previous subsections also apply to the nodal case j = 00 and the cuspidal 
case g2 = g3 = 0. 

Namely, in the nodal case, we get 1-parameter families of algebraically integrable operators 
for r = 1, 2, finite collections operators for r = 4, 5, 7, and conjecturally no solutions for larger r 
(this is confirmed for r < 22). 

In the cuspidal case g2 = 0, (73 = 0, we always have an algebraically integrable operator with 
cyclic symmetry; apart from that, we get 1-parameter families of algebraically integrable ope- 
rators for r = 1,2, and conjecturally no other cases (if r = 3s + 1, this is true by Proposition [32]) • 

4 Operators with several poles 
4.1 Third order operators 

Consider now a third order algebraically integrable operator 

L = d^ + a{z)d + biz) (4.1) 

on an elliptic curve E with several poles zi, . . . ,zn G E. 

It is easy to show that if the gaps of L at a given point are q = r = 1, then the operator must 
be holomorphic at this point, i.e. this case is trivial. So we consider the simplest nontrivial case, 
when the gaps at all the poles are q = r = 2. 

Lemma 4.1. Let L = + a{z)d + b{z) be a Fuchsian differential operator near z = with gaps 
q = r = 2. Let 

a{z) = akz^-"^, biz) = ^ bkz'''^ 

k>0 k>0 

Then L — X has trivial monodromy around for any A if and only if the Laurent coefficients ai, 62 
are zero, and 

bl , , bia^ 

02 = -y, 04 = 04 + —^. 



16 



P. Etingof and E. Rains 



Proof. We have three gaps 2, 2, 4, and three conditions associated to them, which are of degrees 
2, 2, 4. The conditions of degree 2 say that 62 = and 02 = — bf/3. The condition of degree 4 
is hnear in A. The leading coefficient in A is of degree 1, and is a nonzero multiple of oi, so we 
get the condition ai = 0. The constant coefficient is of degree 4 and gives 64 = 04 + B 

Clearly, the same result holds for a Fuchsian differential operator defined near any point 
z = zq. Thus, any algebraically integrable operator (j4.ip with gaps g = r = 2 at all poles would 
necessarily have to be of the form 

/ N \ 3 ^ ^ 

5^ + c - 3^p(z - Zi) a - i^^p'iz - Zi) + ^'ipip{z - Zi), 

\ i=l J i=l i=l 

where pi and c are complex numbers (up to adding a constant). Let us find the conditions on 
the parameters Zj, pi, c for this operator to be algebraically integrable. 

Proposition 4.2. The conditions for algebraic integrability of L are 
c + 3pi = 3^p{zi - Zj), i = l,...,N, 

and 

^{Pi + Pj)p'izi - Zj) = 0, i=l,...,N. 

Proof. The proof is by direct calculation using Lemma |4.1[ ■ 
Corollary 4.3. Let 

^ 3 

F{z,p) = ^pf - - ^{pi +Pj)p{zi - Zj). 

i=l i^j 

Then L is algebraically integrable if and only if {z,p) is a critical point of the function F{z,p) + 
We note that F is the cubic integral for the elliptic Calogero-Moser Hamiltonian 

N 

H2 = ^pf - Yli - ^i)- 

i=l j^i 

Thus, the algebraically integrable operators for a fixed value of c are the critical points of 
H3 + cHi, where Hi = ^pi. 

So, Proposition 14.21 can be viewed as a third order analog of Corollary 12.151 

Remark 4.4. CoroUarv 14.31 is. essentially, a special case of the elliptic analog of Proposition 6 
of P p. 124]. 

4.2 Higher order operators 

We expect that in a similar way one can deal with higher order operators, obtaining families 
of algebraically integrable operators parametrized by critical points of higher Calogero-Moser 
Hamiltonians. Specifically, we expect that if we take the n-th order operator L with indices 
— 1, 1, . . . , n — 2, n and poles zi, . . . , z^v on an elliptic curve, then the algebraically integrable 
operators will correspond to critical points of a degree n elliptic Calogero-Moser Hamiltonian. 
Similarly to Corollary 14.31 this should be a consequence of the methods of [1] and [H] . 
This result obviously has trigonometric and rational counterparts. 

We note, however, that whether the corresponding variety of critical points is nonempty, 
what is its dimension, etc., are, in general, difficult questions. 
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4.3 Operators with symmetries 

It is also interesting to consider operators with symmetries. For example, suppose L is a second 
order operator + u{z) which is even with respect to z. Assume that it has poles at the fixed 
points of 2; — )■ —z (i.e. wq = wi = 1/2, W2 = t/2, and = {1 + t)/2) and at some other 
distinct points izi, . . . , itzTv- Assume that the indices of L at the fixed points Wi are — mj, rrii + l 
for i = 0,1,2,3, and the indices at itzj are —1, 2. In this case, similarly to Proposition 12.15] 
it is easy to show that algebraically integrable operators correspond to critical points of the 
Inozemtsev potential 

3 TV . ^.2 

(see |17t Theorem 0.2] for the case = 1). 

In the same vein, one may consider operators L of order £ = 3, 4, 6 which are invariant under 
the group on an elliptic curve with such Z^-symmetry. Suppose that L has poles at the fixed 
points rjj of Z^, and also at some other points zi, . . . , z^v (taken from distinct Z^-orbits) as well 

as their images under the Z^ action. Let us fix the indices at r]j to be the same as the indices of 

III ■ 

the operator L^ - ^ , where Ij is the order of the stabilizer of and Lqj is a rational homogeneous 



operator (12. ip of order tj with arbitrary integer indices. Also, let us fix the indices at the other 
poles to be —1, 1, — 2, 

Conjecture 4.5. Algebraically integrable operators L as above correspond to critical points of 
the lowest degree {i.e., degree £) Hamiltonian of the classical crystallographic elliptic Calogero- 
Moser system for the group Z^ {with appropriate parameters) defined in [8]. 

Remark 4.6. Conjecture 14.51 may be generalized to the case when the indices of Lqj are not 
assumed to be integers. Namely, in this case we conjecture that operators L with trivial mo- 
nodromy of Lip = Xip around non-fixed points of Z^ (i.e. those for which the monodromy gives 
rise to a representation of a generalized DAHA of rank 1 of type Eq, Ej, Eg defined in [9]) 
correspond to critical points of the lowest degree Hamiltonian of the classical crystallographic 
elliptic Calogero-Moser system for the group Z^ with generic parameters. 

Proposition 4.7. Conjecture 14.51 holds for £ = 3. 

Proof. In the case i = 3, the classical crystallographic elliptic Calogero-Moser Hamiltonian 
of [Hj has the form 

N N 2 2 

H = '^Pi + '^{Arp{zi - r]r)pi + Brp'{zi - T],.)) - '^C^^ p{zi - e'Zj)pi, 
1=1 i=l r=0 i^j s=0 

where t = e := e^'^'/^, p{x) := p{x,t), rjo = 0, 771 = iy/3/3, rj2 = — i\/3/3, and Bi, C are 
parameters. 

On the other hand, consider the Zs-symmetric operator 
2 N 2 



L = d^ + ^{arp{z - 7]r)d + ^rp'{z - r/,)) - 3 X X p(z - e'zi)d 

r=0 i=l s=0 

^ ^ / 3 \ 

+ X X ( - e'^i) + "p{z - e'zi) J . 

i=l s=0 ^ ^ 
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Using Lemma 14.11 and the identity 

p{{l - e^^)z) = — —ip{z) + p{z - rji) + p{z - r/2)), 

we obtain the fohowing conditions for the operator L to be algebraically integrable: 

2 2 
= 3 ^ ^ p{zi - e^Zj) - ^(or - '^)p{zi - r]r), 

jjti s=0 r=0 

2 



^ ^("r - '^)p'iZi - rir)pi + " /3r^ p" (Zi - r]r)^ 



= 3'^'^p'{zi-e''zj){pi + e 'pj), 

for i = 1, . . . , N . But these are exactly the conditions for a critical point of H, with Ar = ar — 1, 
Br = — (3r, and C = 1. ■ 

Remark 4.8. Similarly to the previous subsection, we expect that by considering operators L 
of order ni, n > 1, with symmetry, one can obtain families of algebraically integrable ope- 
rators parametrized by critical points of a Hamiltonian of degree n£ (in momenta) for the 
crystallographic elliptic Calogero-Moser system of [8]. 

It would also be interesting to interpret the complete flow of this system (not only its critical 
points) along the lines of |14j . 



Remark 4.9. Here is a rational version of Conjecture 14.5^ which allows arbitrary i. Namely, 
let 

L = + a2{z)d^-^ + ■ ■ ■ + ae{z) 

be a differential operator with rational coefficients, which is invariant under Z^, such that aj(z) 
vanish at infinity. Let the nonzero poles of L be zi, . . . ,zn (taken from different Z^-orbits) as 
well as their images under the symmetry. Suppose that L has arbitrary integer indices at 0, and 
indices —1, 1, ...,£ — 2,£ at zi, ... , zn- 

Conjecture 4.10. Algebraically integrable operators L with such properties correspond to critical 
points of the rational Calogero-Moser Hamiltonian {of degree i) for the complex reflection group 
Sn X Zf (see e.g. [8]). 

It follows from the above that this conjecture holds for i < 3. 

Moreover, we expect that a similar conjecture holds for operators of order ni with Z^- 
symmetry. Namely, in this case we should require that the indices at are those of Lq, where 
Lq is a rational homogeneous operator of order i with integer indices, and we conjecture that 
algebraically integrable operators are parametrized by critical points of the higher order rational 
Calogero-Moser Hamiltonian (of order ni). 
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